Decay of correlations in the dissipative two-state system 
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We study the equilibrium correlation function of the polaron-dressed tunnelling operator in the 
dissipative two-state system and compare the asymptoptic dynamics with that of the position corre- 
lations. For an Ohmic spectral density with the damping strength K = |, the correlation functions 
are obtained in analytic form for all times at any T and any bias. For K < 1, the asymptotic 
dynamics is found by using a diagrammatic approach within a Coulomb gas representation. At 
T — 0, the tunnelling or coherence correlations drop as t~ 2K , whereas the position correlations show 
universal decay oc t~ 2 . The former decay law is a signature of unscreened attractive charge-charge 
interactions, while the latter is due to unscreened dipole-dipole interactions. 
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The simplest model that allows to study the interplay of tunnelling and dissipation is the spin-boson model . 
Despite its simplicity, it exhibits generic features of many complex systems in physics and chemistry and has found 
widespread applications. It has been adopted to describe diverse systems like the tunnelling of atoms in atomic- force 
devices ||, or the dynamics of the magnetic flux in a rf-SQUID Q, just to mention a few. 

The Ohmic spin-boson model exhibits a dynamical phase transition between coherent and incoherent tunnelling. 
Both for the expectation value (a z (t)) (the population) and the a z autocorrelation function, the transition from 
oscillatory to overdamped behaviour was found to occur at a damping strength K = | at T = and zero bias 
Concerning the dynamics at long times, these functions behave quite differently. The factorized system-reservoir initial 
state for (a z (t)) leads to exponential decay J^H, whereas the correlated initial state of the symmetrized equilibrium 
correlation function implies an algebraic decay oc 1 /t 2 |8|,|| . 

The dissipative two-state system (TSS) is described by the spin-boson Hamiltonian 
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where Ho = —U (Aa x + ea z ) /2 describes the isolated TSS with level splitting TiA and a bias energy he. The eigenstates 
of a z are the two localized states \R) and \L) at positions ±a/2. The effects of an Ohmic bath are captured by the 
spectral density J(u>) = (2irhK/a 2 ) uj e~"/" c , where K is a dimensionless damping strength and uj c is a cutoff for 
the bath modes. Here we are interested in the scaling limit A t /uj c — > 0, where A r = A(A/u> c ) K ^ 1 ~ K ^ is kept finite. 
In this limit, the population (a z (t)) and the a z autocorrelation function are universal functions of the renormalized 
frequency A r , i.e., there is no explicit dependence on oj c . 

Recently, there has been considerable interest for observables of the tunnelling operator <r x |l^,[ll|] . The correlations 
of cr x describe the response of the system to a change in the barrier width or in the barrier height. In atomic force 
devices, this can be contrived by modulating the distance between surface and tip, and in the rf-SQUID by varying 
the critical current [Q. It was shown that, because of an explicit ui c dependence, both the expectation value and the 
autocorrelation function of the bare a x vanish in the scaling limit . 

In this Letter, we study the correlation function of the polaron-dressed tunnelling operator a x = XJa x XJ~ x which has 
a nontrivial scaling limit. The polaron transformation U — exp{— ia z ^2 a s a p a /2h} induces adiabatic displacements 
{s a } = {aCa/rriaUJa} of the bath modes during the tunnelling process (cf. Ref. [jl|l3|). The dressed tunnelling 
operator, a x = |i?)(L| exp{— i ^2 Q s a p a /h} + h. c. , transfers the particle from one localized state to the other and 
simultaneously shifts the bath modes by the displacements {±s Q }. It is just in this way that coherent oscillations are 
induced and therefore a x is called the coherence operator. 

The quantities of our interest are the equilibrium correlation functions of the position, C z (t) = (c z (t)a z (0))p— {o z ) 2 p, 
and of the coherence operator, C x (t) — (^(^^(O))^. Our emphasis is put on their asymptotic behaviour for general 
K < 1. Using a real-time path integral approach, an exact formal series expression for the symmetrized position 
correlation function S z (t) = Re C z (t) has been derived in Ref. [Fbil , 
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Similarly, the respective response function XzW — — ®K) I m C z (t) i s found as 



„ OO OO „t 

Xz(t) = --EE (~A 2 )' i+m tan(TrA) / P&U^te} ]T & & G «,™^ ■ (5) 

Here we have put A 2 = A 2 cos(irK)/2. The integration symbol contains the time-ordered integrations over k flips in 
the negative and I flips in the positive time branch, 
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In these expressions, the summations over the intermediate diagonal (sojourn) states have already been performed. 
Further, the £j summation takes into account the two possibilities for each off-diagonal (blip) state the system takes. 

The effect of the bias is in the factors £>i + m = cos(e Y^j=" 1 £j T j) an d Dn.m — sin(e X}j=i™ £j T j)i where Tj — t%j — ^2j-i 
is the length of blip j. The bath-induced correlations are expressed in terms of the interaction S Ptq = S(t p — t q ), where 
S(t) = 2A ln[ (h(3uj c /n) sxnh^t/hfl) ]. In the equivalent Coulomb gas picture, each factor exp[— S(t p — t q )] represents 
the interaction of a neutral pair of charges, which we refer to as a dipole. The full intra- and interblip correlations 
are in the factor 
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For C x (t), the procedure is more complicated, since the polaron-dressed operator a x acts on the full system-plus- 
reservoir space. Therefore, we have reconsidered the elimination of the bath modes. We have found that the resulting 
influence functional can be cast into the standard Feynman- Vernon form at the expense of introducing modified system 
paths. The findings can be put in simple terms: the system's jumps at times zero and t enforced by the operation of 
a x do not give rise to bath correlations in the influence functional. 

The contributions to C x {t) can be divided into two groups. Either the system jumps from a sojourn to a blip state 
at time zero (group A) or vice versa (group B). This extra jump is accounted for by the upper index r — 1 in the 
integration element (||). For both groups, the system is finally in a blip state. The symmetrized correlation function 
5a; (t) = KeC x (t) and the response function Xx(t) = —(2/h) Q(t)lmC x (t) are found as 
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The superscripts {. . . } A and {. . . } B indicate that the blip labels are subject to the constraints 

in+m = £n+i (group A) , £ n+m = -£„ (group B) . (13) 
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The interaction factors G^ m and Gf t m differ from the form (Q) for G„ 
and t' — t. We have 



by the removal of the charges at times t' = 
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where A^ k and A^ fe describe the interblip correlations for the modified sequences of charges. If j, k ^ n + 1 and 
7^ n + to, A^ fc is again given by (||). In all other cases, the interactions of the missing charges in Eq. (||) have to be 
dropped. For instance, for j = n + 1, we have A^ +1 k = 52n+2,2fc-i — S l 2n+2,2/c- There are analogous modifications in 
A^ fc for j, k = n and n + m. For instance, we have A^ fe = S 2n -i,2k - S 2n -i,2k-i- 

For the value K = i , the above series can be summed in analytic form using the concept of collapsed dipoles jl4| . 
Putting K — i — k with k < 1, the phase factor cos(irK) w ttk vanishes in the limit k — > 0. In order to have a finite 
contribution for K = |, each factor cos^if) has to be compensated by the 1/k "short-distance" singularity arising 
from the breathing mode integral of a dipole with interaction e~ s ^ ~ {uj c t)^^ 2k \ The combined expression is 
termed a collapsed dipole, yielding the finite contribution 
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A collapsed dipole does not interact with other charges and it is insensitive to a symmetric bias factor. In contrast, an 
odd bias factor in Eq. (16) prevents a dipole from collapsing. However, to give a nonzero contribution, the extended 
dipole must be free of a cos(irK) factor. Within an extended blip [sojourn] state of length r [s], the system may make 
any number of visits of duration zero to a sojourn [blip] state. This is described by the insertion of a grand-canonical 
gas of non- interacting collapsed sojourns (CS) [or blips (CB)], yielding a factor e _7r / 2 [or e _7S ]. Note that there is 
a multiplicity factor 2 for collapsed blips. As a general rule, an extended sojourn interval, say is free of insertions 
only if the subsequent blip is weighted with a factor £fc+i. In the expression (|j|), e.g., the first sojourns in each time 
branch are without insertions. Employing these concepts, the following results for the various correlation functions 
are obtained. The a z correlation functions are given by Q 



S z (t) = e-* - F 2 + {t) - F 2 (t) , *.(*) = m)F+(t) e^ t/2 , 
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with f+{x) — cos(ir), f-(x) = sin(x). A diagrammatic representation of S z (t) and Xz{t) is displayed in Fig. |]. 
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FIG. 1. The diagrams for \z{t) (right) and the contribution Q s {t) to S z (t) (left). The full and dashed lines represent sojourns 
and blips, respectively. An empty box represents the insertion of a CS gas within a blip interval, the full box stands for a CB 
gas inside a sojourn interval. The upward and downward spikes symbolize charges. 



In the asymptotic regime T = and t 3> I/7, we find from Eq. (|l7| ) 

/4\ 2 / 7 2 \ 2 1 , x 16 7 2 cos(rt)e- 7t / 2 

The exponential decay of the response function is due to the grand-canonical sums of collapsed dipoles in each interval 
(except for the first sojourn). In contrast, the function S z (t) decays algebraically because the first sojourn in the 
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positive time branch is free of collapsed blips. Therefore, this interval gets very large and is effectively limited by the 
overall length t. The l/t 2 law reflects the bare interaction between the two dipoles displayed in Fig. [j] (left). 

Now switch to the a x correlation functions. Since there is one cos(ttK) factor more than dipoles in Eq. @, S%(t) 
vanishes as K — > 4. Therefore, only group A contributes to the symmetrized correlation function, yielding the damped 
oscillatory behaviour 

S x (t) = cos(ei) e~ 7 * /2 . (20) 
The corresponding diagram is sketched in Fig. 0. 



. £j . 

-oo t 

FIG. 2. The diagram describing S x (t). The bullets mark transitions which are free of bath correlations because of the 
modified influence functional. 
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FIG. 3. The contribution of group A (left) and group B (right) to Xx(t)- Each diagram has only one extended dipole. 

For the response function, the contributions of both groups are shown in Fig. S The expressions are combined to 

Xx {t) = (2/ft)[cos(ei)F+(t) + sm(et)F_(t)] . (21) 
At T = and times t > I/7, we find from Eq. (j|l]) 
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The algebraic decay law arises from the contribution of group B. According to the above rule, the sojourn interval in 
the positive time branch is free of a CB gas factor and therefore is effectively very large for t 3> I/7. The l/t law in 

S W rv 1 /* for K = i. 



Eq. ( |22| ) is simply the signature of the bare intra-dipole interaction, e~ ' ' oc l/t for K 

For general K, it is not possible to sum the series for the a z and a x correlation functions in analytic form. 
Nevertheless, it is possible to extract the long-time behaviour of the correlations. The essential modifications concern 
the CB and CS gas factors inserted in a given interval. For K ^ h, the dipoles are no longer collapsed and thus the 
respective grand-canonical sum can not be performed in analytic form. However, because of the alternating sum, the 
charges which are partitioned off by the weight factors £k+i form clusters of effective length 1/A r . 

Consider first the a z correlation functions at asymptotic times t 3> 1/A r . For x z (i), there is a single neutral cluster 
surrounding the origin of the time axis. Therefore Xz(t) decays exponentially. In S z (t), we have a neutral cluster in 
each time branch. Since in both branches the initial sojourn is free of insertions (see above rule), the two clusters 
are near the origin and near t, respectively. The correlations between the two neutral clusters are dominated by 
the bare dipole-dipole interaction which is —2K/t 2 at T = 0. The interaction has a universal power-law form which 
is independent of the coupling strength. The power 2 is a characteristic feature of Ohmic dissipation. The dipole 
moment emerges to coincide with the static susceptibility S . In the end, we have 
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S z {t) = -2K[h X ^/2] 2 ^. (23) 
In the frequency domain, Eq. (|23|) corresponds to the generalized Shiba relation 0] for the spectral function, 
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Next, consider the a x correlation function for t 3> 1/A r . In group A, there is only a single neutral cluster near 
to the origin of the time axis. Hence, both S x (t) and Xx (*) decay exponentially at asymptotic times. In group B, 
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we have a charged cluster in each time branch, satisfying overall neutrality. Again, as a result of the above rule, the 
clusters are situated near the origin and near t and therefore are roughly separated by an interval of length t. The 
clusters interact with unscreened charge-charge interaction e~ s ^ oc t~ 2K . This interaction directly determines the 
long-time behaviour of it) and (t) ■ The contributions of group B predominate over the exponential contributions 
of group A for t 3> 1/A r . Thus we find the asymptotic behaviour 

S x (t) oc e- s ^ oc t~ 2K , Xx {t) oc e" s « oc r 2K . (24) 

We note that the prefactor in the algebraic decay law for S x (t) vanishes accidentally as K — > i. The results (|2~i| ) 
have been indicated numerically in Ref . . The slow decay of \x (t) for K < i implies that the static susceptibility 
diverges algebraically, oc T 2K ~ 1 . Interestingly, this regime coincides with the coherence regime for the population 
(a z (t)) at zero bias p|,[TT[| . 

In conclusion, we have given exact formal expressions for the tunnelling or coherence correlations in the dissipative 
two-state system for general damping strength K and we have presented results in analytic form for the particular case 
K = i . The differences in the long-time behaviours of the position and coherence correlations have been illustrated 
in terms of a simple charge picture. 
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